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Abstract. We study the difference between verticai and conicai square functions 
in the abstract and also in the specific case where the square functions come from 
an elliptic operator. 



o 

. 1. Introduction 



The purpose of this article is to draw attention to differences between vertical and 
Q ■ conical square functions. By vertical square functions, we mean the usual Littlewood- 

! Paley-Stein functionals. By conical square functions, we mean the area functionals of 

Lusin type. Our interest in this subject was triggered by the recent unpublished work 
of Dragicevic and Volberg jDV] . Let us first describe what they proved. 
<^ Let A = A(x) be an n x n matrix of complex, L°° coefficients, defined on M n , and 

! satisfying the 

jh ; m 2 < Re A£ • £ and \A£ • C| < A|£||C|, 

for £, £ G C n and for some A, A such that < A < A < oo. We define a second order 
divergence form operator 

L/ = -div(AV/), 

which we interpret in the sense of maximal accretive operators via a sesquilinear form. 

Proposition 1.1 ( [D Vj ) . If A is real and 1 < p < oo, there is a dimension free 
bilinear estimate 



OO 



1.1) // _ |V B e-* L /(y)||V y e-^(y)|dydt<C(p,A,A)||/|| p ||p||, 



n+l 



Here p' is the conjugate exponent to p. 

Set aside the dimension free bound, this result is striking in view of the following 
vertical square function estimate. 

■ Proposition 1.2 ( |Ausj ). If A is real and 1 < p < q+(L), 

r / roo \ p/2 

(1.2) / / \V y e- tL f(y)\ 2 dt) dy<C(p,n,X,Ay\\f\\l. 

Furthermore, this estimate fails for p > q+(L) (if q+(L) < oo). 
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The number q+{L) has been introduced in [Aus] . as well as the three other numbers 
p_ (L) , g_ (L) and p+(L) as limits of the following intervals. The interval (p_ (L),p + (L)) 
is the maximal open interval where the heat semigroup {e~ tL }t>o is uniformly bounded 
on L p or equivalently the semigroup satisfies LP — L q off-diagonal estimates when 
P-(L) <p < q < p+(L) , see (13. ip below. Analogously, (q_(L),q + (L)) is the maximal 
open interval where {\/tVe _<L } t> o is uniformly bounded on LP or satisfies LP — LP off- 
diagonal estimates when q~(L) < p < q < q+(L) . These intervals also determine up 
to endpoints the range of LP boundedness of the functional calculus, Riesz transform 
and vertical square functions. One has p~(L) = q~(L), (q + (L))* < p+(L) — where 
q* = qn/(n — q) when q < n and q* = oo otherwise. Also, P-(L) = q~(L) = 1 and 
p+{L) = q+(L) = oo if n = 1. For n = 2, or for n > 3 and L with real coefficients, the 
same is true except that one can only say that q+{L) > 2 for n > 2 and this is sharp. 
Additionally, if n > 3 and A with complex coefficients, then p~(L) < 2n/{n + 2) 
and p+(L) > 2n/(n — 2). See |Aus] for full details. Also |HMMc] proves the latter 
inequalities to be sharp using an example of Freshe |Frej . 

Since q+(L) can be arbitrary close to 2, one cannot deduce the bilinear estimate in 
Proposition 11.11 from the vertical square function estimate. So the bilinear estimate 
seems to exhibit some special feature that the vertical square function does not have. 
Indeed, bilinear integrals as above can also be estimated using conical square functions 
thanks to an averaging trick that appears in |FS] and |CMSj : 

F M G M dydt = K > I (// r „ +1 FMG M h(^) ix 

with h the indicator function of the unit ball and b n its volume, so that 

< Cn \\ShF L \\ShG L 



II F(y,t)G(y,t)dydt 

J JR1+ 1 



with 



S h F(x)=[ I \F(y,t)\ 

'\x-y\<y/i 



2 dydt 



1/2 



Hence, applying this to F(y,t) = V y e~ tL f\y) and G(y,t) = V y e~ tL g{y) it becomes 
natural to expect the corresponding conical square function estimate holds in a larger 
range of p than the one for (II. 2p . Indeed, we shall show as part of Theorem 13. II 

Proposition 1.3. Lf A is real and 1 < p < oo, 

(1.3) / (If \V y e- tL f(y)\^Y 2 dx<C(p,n,X,Ar\\f\\ P P - 

JR" \JJ\x-y\<Vi 1 ' J 

Thus ( II. ip holds at least with a dimension dependent bound . We shall also study 
(II. 3p for all complex A and show it holds when p~(L) < p < oo and fails when 
p < p~{L) (if P-{L) > 1). This is consistent as p~(L) = 1 when A is real. This also 
improves |Aus[ Corollary 6.10] where ( II. 3p was obtained in the range p~{L) < p < 
q+{L). The bilinear inequality as in Proposition 11.11 then holds for a restricted range 
P-(L) <p <P-(L)' = p+(L*). 

This leads us to the main point of this article about comparison between vertical and 
conical square functions. Propositions 11.11 and 11.31 show that the ranges of p below 
2 are the same but differ above 2. One may wonder whether there is an abstract 
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principle behind this. But this is not the case. Aside from p = 2 for which the 
averaging trick yields that they are equivalent, vertical and conical square functions 
only compare one way for p ^ 2 in the sense that one is automatically controlled by 
the other and simple examples show the converses fail. More precisely, for p > 2, it 
is well-known and we shall recall why in Section [21 that a vertical square function 
controls the corresponding conical one. We shall also prove, and it seems this is not in 
the literature, that for p < 2, the conical square function controls the corresponding 
vertical one. Comparing the ranges for (11.31) and (II. 2p already furnishes a counter- 
example for the converse in the p > 2 range and an example where the converse holds 
in the p < 2 range. We note this can be done on a space of homogeneous type. We 
shall also study some weigthed comparisons using extrapolation. 

We finish this introduction by the following observations. As explained before the 
range of p for (ll.2p is tight to the range of LP boundedness for ^/tVe~ tL . As the 
p < 2 range for (jl.3p is a priori smaller than or equal to the p < 2 range for (II. 2p , we 
obtain the best possible result by showing they are equal. For p > 2 we exhibit a new 
phenomenon. 

Our results show that the p > 2 range for (II .3p is linked to the rate of decay in 
the L 2 off-diagonal estimates. If the latter is fast enough then one obtains the full 
range (2, oo) as it is the case in (II. 3p . In fact this LP estimate amounts to proving 
boundedness of some vector-valued operator from LP into the parabolic version of 
the tent space T% of Coifman, Meyer, Stein |CMSj . When the L 2 off-diagonal decay 
is fast enough we can prove, basically following the Fefferman- Stein argument, that 
\V y e~ tL f(y)\ 2 dydt is a parabolic Carleson measure for / G L°°, which is nothing 
but an L°° — > T£° estimate. One can then interpolate for 2 < p < oo. When the 
rate of decay is slow (for example polynomial with small exponent) this argument 
does not seem to adapt and one needs other tools. This is the case for the conical 
square function based on (p(t 2 L) when A is complex and (p not smooth as the origin. 

1/2 

An example is the Poisson semigroup since (p(z) = e~ z in this case. A different 
ingredient then comes into play, which is the decay at of ip or the order at which 
it vanishes, combined with the definition of p + (L). For instance, in Section [3T41 below 
we shall prove the following and this is the p > 2 range that is interesting for our 
discussion here. 



Proposition 1.4. For m a non negative integer and f G LP , then 

-tL 1 / 2 

\x-y\<t 



,dydt\ p ^ 



(1-4) / // \tV v 4(1?L) m e- tL 'f){y)\ 2 -^) dx<C\\f\\l 

whenever 

(1.5) P~(L) < p < 



np + (L) 



n 



{2m+l)p+{L)' 



For (2m + l)p + (L) > n then the right hand side in (II. 5p becomes oo. 

If L had been the Laplacian, 2m + 1 would just be the number of vanishing mo- 
ments for the kernel of the convolution operator V y ,t(t 2 L) m e~ tLl/2 or the number of 
derivatives in front of the semigroup. Here we have 2mth order "vanishing" coming 
from the exponent of the second order operator L, and +1 comes from the gradient. 
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A few comments are in order. We first point out that, since vertical and conical 
square functions are equivalent in L 2 , one may integrate by parts in t and use prop- 
erties of the semigroup to pass from any choice of non-negative integer m to another, 
in the case p = 2. In fact, one may even take m of the form m = k + 1/2, with k a 
non-negative integer. For p > 2, different values of m apparently need no longer be 
equivalent; the conclusion of the proposition yields a better range of p for larger m 
(up to the critical value with (2m + l)p + (L) = n). 

In particular, the case m = of (II. 4p gives standard area integral estimates for 
weak solutions of the equation 

(1.6) d 2 t u + div AVu = 0. 

When A is real (in which case P-(L) = 1, p+{L) = oo)), such estimates may be 
obtained as a special case (the "block matrix" case), of the result of Dahlberg, Jerison, 
Kenig in jDJKj . using the fact that one has non-tangential estimates for the solutions 
in every L p , 1 < p < oo. The present argument allows for a direct 
(and simpler) proof than that in [DJK] in this special case. Moreover, it has the added 
virtue of applying to the case of complex coefficients. Of course, we do not address 
the question of "full" coefficient matrices (i.e., those that need not be in block form), 
as is done in [DJK] . 

Acknowledgments. This work was started years ago while the authors were all 
visiting the Universidad Autonoma in Madrid on the occasion of a special program in 
harmonic analysis and, after a latency period, finished this year while the authors were 
all visiting the Center for Mathematics and Applications of the Australian National 
University. We are very grateful to these institutions for hospitality and financial 
support. Also we want to express our thanks to O. Dragicevic and A. Volberg who 
showed us their unpublished work. 

2. Vertical versus conical 
For a locally square integrable function / on R" +1 , denote 

sm=(ff l/(y,t)| 2 ^) 1/2 , xer. 

and 

Vf(X)= (/>o' /(M)|2 t) ' XeRr 
Sf is lower semi-continuous hence a measurable function. Measurability on V f follows 
from the local square integrability of /. 
We remark that 

\\Sf\\l = b n \\Vf\\l 
with b n the volume of the unit Euclidean ball. 

2.1. Comparison in Lebesgue spaces. 

Proposition 2.1. Let f be locally square integrable on 

(a) For 2 < p < oo , 

\\Sf\\ p <C(p,n)\\Vf\\ p . 
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(b) For0<p<2, 

\\Vf\\ p <C(p,n)\\Sf\\ p . 

(c) The converses fail for all p ^ 2. 

Proof. Part (a) is standard and appears already in (StJ p. 91]. For the sake of self- 
containment, we recall the argument. As p > 2, q = p/2 > 1 and we can estimate 
ll^/llp = II (Sf) 2 \\p/2 by dualizing against a function h £ L q . Now, the averaging trick 
applies and yields 



(SfY(x)h( X )dx = b n \f(y, t )\ 2 [-—— / h(x)dx 



(y,t) 



dydt 



' \\B(y,t)\ Jb 

<b n [ (Vff(y)Mh(y)dy 
< b n \\V ff p \\Mh\\ ql 

and one concludes using the boundedness of the maximal operator over balls M in 
L q ' . 

Let us now prove Part (b). Fix < p < 2, / with Sf £ L p and A > 0. Since 
x i — y Sf(x) is lower semi-continuous, the set O = {Sf > A} is open. Let F be the 
complement of O in M. n , R(F) be the union of the cones \x — y\ < t with vertices 
x £ F. We also set 6 = {x £ R : M(xo)O) > 1/2} and F = M n \ 6. We 
note that O C O since O is open, and thus F d F . If y £ F and t > we have 
|OD B(y,t) |/|£(y,t)| < 1/2 and consequently, \F n B(y,t)\/t n > 6„/2. Hence, 

2 dxdydt 
t n+1 



(Sf) 2 (x)dx= I // 

xGf </ </|z-j/|<t 



FnB(y,t)\ t ., 2 dydt 



R(F) t 

b n f f „ j.m2 dydt 



> f / . / l/(v,t)l 



2 
6, 



</t>0 



(V/) 2 (|/)dy 



>yA 2 |W>A}nF|. 



Besides, for < r < p, using the weak type (1,1) for M, 

~ 2 ■ S n f 

\{Vf > A} n 0\ < \0\ < 2 ■ 3" \o\ < — J o (SfY(x) dx. 



Hence, 



f* POO 

/ (Vfy(x)dx=p / X?- 1 \{Vf> \}\d\ 

JR n JO 



2p 



oo 



< ^ / A^ 2 " 1 / (SfY(x)dxdX 

"n Jo JSf<\ 

/*oo p 

+ 2p-T X^- 1 (Sf) r (x)dxdX 
Jo Jsf>\ 
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2p_ 2p^\ r {Sfn )d 

We now finish the proof with Part (c). It is convenient to introduce 

Sf(x)=[[ \f{v,t)\%£, Vf(x)=[ 

J J\x-y\<t 1 Jt>0 1 

Note that Sf = S(\f\ 2 )^ 2 and Vf = V^(|/| 2 ) 1/2 , so that, for a locally integrable 
function / on W r l +1 , we seek to disprove the inequalities 

(2.1) \\Sf\\ p <C\\Vf\\ p , 0<p<l; 

(2.2) \\Vf\\ p <C\\Sf\\ p , Kp<oo. 

We write the argument so that it is easy to adapt it to a space of homogeneous type, 
denoting v(B) the volume of a ball and using implicitly the doubling property in the 
argument. See Remark 12.21 below. 

For ( 12.1 p we consider f N (x,t) = A r_1 1 x( x ) Xo(t/N) with N ^> 1 and where x is 
the characteristic function of the unit ball -8(0, 1) and Xo denotes the characteristic 
function of the interval [0, 1]. On the one hand, 

Vf N (x)= / \f N ( x ,t)\- = N- 1 X (x) / dt = X (x) 
Jo t Jo 

and therefore 

||V7*||5 = v(B(0,l)). 

On the other hand, fixed |x| < N/8, if \y\ < 1 we have \x — y\ < N/4 (provided N > 8) 
and then 

Sf N (x) = N~ l [[ tx(y)Xo(t/N) 



= N~' 

> N" 1 

> C 



\x- y \<t tv(B(y,t)) 
dtdy 



\y\<\ J\x-y\<t<N v(B[X, t)) 

dtdy 

N/A<t<N v(B{X, t)) 



v(B(0,l)) 



This implies 



v(B(0,N)Y 

v(B(0,l)Y 



\\Sf N \\ p p > [ Sf{xYdx>C 

J\x\<N/8 



I\x\<n/8 V 7 " v(B(0,N))p-i' 
Gathering the obtained estimates 

11^16 >c v(B{0,N)) 1 -* 

\\vf N \f p - v(B(o,i)y-p- 

Thus (12. ip cannot hold as v(B(0, N)) increases to oo and 1 — p > 0. 

For (12. 2p we consider f^(x,t) = t v (B(x, t)) Xn(x) Xo(t) with N ^ 1 where xn is 
the characteristic function of the ball of radius 1/N. We first calculate V /n- 

r°° dt 

vf N (x)= / \f N (x,t)r 

./n t 
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> Xn{x) [ v(B(x,t))dt 

J 1/2 

>C X n(x)v(B(x,1)) 
>C X n(x)v(B(0,1)) 



and therefore 



\\Vf N \\l>Cv(B(0,l/N))v(B(0,l)Y. 

We find an upper bound for ||S/jv||p- We notice that if \x\ > 2, \y\ < 1/N and 
< t < 1 we have \x - y\ > 1 > t (if N > 1). Thus, Sf N (x) = if \x\ > 2. On the 
other hand, for all x G lR n : 

Sf N (x)= [[ X N(y)dydt<v(B(0,l/N)). 

J J\x-y\<t<l 



Then, we obtain 



so that 



\\Sf N \\ p p = [ Sf N (x) p dx< Cv(B(0, l/N)Yv(B(0,l)) 

J\x\<2 



\\sf N \\i <c v(B(o,i/N)y-i 



\\Vf N \\ p p - v(B(0,1))p-i 
which goes to as N — > oo if p > 1. □ 

Remark 2.2. The reader can notice that this theorem generalizes to spaces of homoge- 
neous type X with infinite volume and at least one point that is not an atom (which 
plays the role of 0). That is is changed to X x R + and in the definition of Sf(x) 



the measure has to change to 



dfi(y)dt 



2.2. Weighted estimates via extrapolation. Let us present a weighted version of 
Proposition [2TTJ using extrapolation. That is, L 2 estimate with suitable Muckenhoupt 
weights imply L p comparisons in weighted spaces. Let A p , 1 < p < oo, denote 
the classical Muckenhoupt classes of weights and RH p , 1 < p < oo, the class of 
reverse Holder weights. See for example |AMj . Again everything extends to a space 
of homogeneous type as in the remark above. We stick to the Euclidean space for 
simplicity. 



Proposition 2.3. Let f be a locally square integrable function on 

(a) For 2 < p < oo and w 6 A p / 2 

\\Sf\\ LP(w) <C(p,w) \\Vf\\ LP(w) . 

(b) For < p < 2 and w G RH^/ P y 

\\Vf\\ LP{w) <C(p,w) \\Sf\\ LP(w) . 
Proof. We begin with Part (a). Given any w G A ra we easily have 

(2.3) \\Sf\\l Hw) = [ [[ \f(y,t)\^w(x)dx 



pn+1 



\x-y\<t 

\fM\ 



2 w(B(y, t)) dydt 
\B(y,t)\ t~' 
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We note that if w G Ai, that is, Mw(y) < w (y) f° r a - e - 2/ £ K"; then we have 
for all £ > 

w(B(y, t)) ^ ^ ^ a e ^ ^ R „ 



Then, we obtain 



\\Sf\\h M <b n [w] Al / / \f(y,t)\\w( y )dy = b n [w] Al \\Vf\\% {w) . 

JR n JO 1 

Next we invoke the Rubio de Francia extrapolation theorem (see [Rub] . |Garj for 
the original result, and |CMPlj . |CMP2j for a statement written in terms of pairs 
of functions) for the pairs ((S f) 2 , (V f) 2 y. the starting estimate in L^iw) for every 
w G A\ implies that for every 2 < p < oo and w G A p /2 

\\Sf\\ LPH <C{p,w) \\Vf\\ LPH . 

Strictly speaking, the argument applies whenever the left hand side is finite. This is 
the case if / is a priori bounded with compact support in Monotone convergence 

implies that the inequality is valid for all locally square integrable function /. 

For the reverse estimate in Part (b), we recall that w G RH^ if for every ball B we 
have 

w(x) < [w} RHoc — / w(y) dy, a.e. x G B. 



\B\ JB 

Then, using Lebesgue's differentiation theorem we obtain that for a.e y G W 1 and for 
all t > 

w(y) < sup — - — - I w(x) dx 

0<r<t \B{y,T)\ J B (y,r) 
< Hflffoo SU P in/ \ I / id/ TTi / W{z)dzdx 

o<T<t\B[y,T)\ Jb( v ,t) \B[y,t)\ J B (y,t) 

~ [w]rh ~ "W 

Thus, for every w G RHoo by (12. 3 p we have 

,dt 



\\Vf\\b W = /__ / \f(y,t)\ 2 jw(y)dy 



w(B(y, t)) dydt 
\B(y,t)\ r 



/* poo 

< / / l/(v,f)l 5 

JK" Jo 
= WRH^b- 1 \\Sf\\ 2 L 2 iw y 

Considering the pairs (F, G) = ((Vf) 2 , (Sf) 2 ) we have obtained that 

/ F(x) w(x) dx < [w] RHoo b~ l / G(x)w(x)dx, V 'w G RH^. 

JR n JR n 

We take an arbitrary p with < Pq < 1 and set go = r — 1- Then the last estimate 
holds in particular for every w G A r / po fl RH/ qo / r y. We apply the extrapolation 
theorem for limited ranges |AM[ Theorem 4.9] (see also |CMP2j ) to conclude that 
for all po < q < qo 



/ F(x) q w(x) dx < C w I G(x) q w(x) dx, Vw G A q/po n RH 

JR n JR n 
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whenever the left hand side is finite. This is the case when / is bounded with compact 
support in and can be removed by monotone convergence to allow all locally 

square integrable function /. Next, we fix < q < 1 and w G RH(i/ q y. Then, 
w G Aoo and there exists < po < q such that w G A q / P0 . Thus we can apply the last 
estimate since 0<po<q<l = q and w G A q / po D RH( qo / g y. Hence we have proved 
that for every < p < 2 and w G RH( 2 / P y 

\\Vf\\ LP{w) <C(p,w) \\Sf\\ LP{w) . 

□ 

Notice that from the argument one sees that the extrapolations take initial estimates 
in L l (w). Indeed, from the beginning one could have worked with the operators S and 
V defined above. The argument just presented shows that for every locally integrable 
function / on R" +1 , if 1 < p < oo and w G A p then 

\\Sf\\ LPH <C{p,w) \\Vf\\ LP{w) , 

and if < p < 1 and w G RH^/ p y, then 

\\Vf\\ LP{w) <C(p,w) \\Sf\\ LP{w) . 

3. Square functions for typical functions of L 

Consider the operator L defined in the Introduction. We introduce the following 
conical and vertical square functions 

1 /2 

GpU)(v) = (j[ >o l«V B , ( e- li "V(!/)rf ) , 

G h (f)(y)=( I \V y e- tL f(y)\ 2 dt 
\Jt>o 

The P subscript refers to the fact that we are dealing with the Poisson semigroup 
e~* Ll/2 for L. The h subscript refers to the heat semigroup e~ tL . The curly letters 
are for the conical square functions and the capital letters for the vertical ones. So 
from our general observations we know that ||£7p(/)|| p < II^p(/)|Ip for 2 < p < oo and 
||Gp(/)||p < ||^p(/)||p for < p < 2 and similarly for the heat versions by making a 
change of variables in t. Note that these square functions all contain a spatial gradient. 
Hence we are not working within the functional calculus of L. 

We want to compare the LP norms of each square functions with the LP norm of 
the original function /. 

For p = 2, a mere integration by parts (see |Aus] ) yields that 

l|Gp(/)|| 3 +||C? fc (/)|| a «C7(A,A)||/|| a . 
As seen above, conical square functions behave as the vertical ones in L 2 . 
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We turn to a summary of results on LP. Let p* = if p < n and oo otherwise. Let 

j r n—p r 

us remind the reader that the exponents p±(L), q±(L) were defined in the introduction, 
in the discussion following the statement of Proposition 11.21 

Theorem 3.1. (1) Gh is bounded on L p for p~(L) < p < q + (L). 

(2) §h is bounded on LP for p~(L) < p < oo. 

(3) Gp is bounded on LP for p~(L) < p < q+(L). 

(4) Qp is bounded on L p for p_(L) < p < p + (L)* . 

The upper bounds are optimal except maybe for Gp- The lower bounds are all opti- 
mal. 

The converse estimates \\f\\ p < ||^(/)||p are valid for all p £ (l,oo) and f £ L p r\L 2 
and all four square functions. Hence, each defines a new norm on L p for p in the 
corresponding range above. 

Fix /i £ (0,7r/2) and 1 < p < oo. We say, following [Ausj . that a family of linear 
operators (T z ) ze ^ satisfies L p — L q off-diagonal estimates if there exist constants 
c,C such that for all z £ E M := {z £ C*; | argz| < /i}, all Borel sets E, F and all 
/ £ L P (E), we have 

cd(E,F) 2 

n r 1 1 \ - ' - 

(3.1) \\T z {f Xe)\\l^f) <C\z\-^- P -^e 1*1 \\f\\ LP{E) . 

This holds for fi < n/2 — uj with u the type of L, T z = (zL) m e~ zL with p~(L) < 
P < q < P+(L) and T z = \z\ l/2 W(zL) m e~ zL with q-(L) = pJ(L) < p < q < q + (L), for 
any non-negative integer m. See |Aus| Chapter 3]. 

For 1 < p < oo, we recall that the tent space Tf denotes the space of locally square 
integrable functions in such that Sf £ L p (W l ) with the notation of Section [2] 
The norm in T^ is given by US/Hp as defined in Section |2] Note that changing the 
aperture of cones yields equivalent norms. For p = oo, we let T 2 °° be the space of 
locally square integrable functions in such that 

n/ik~ = (4 // i/(y,t)i 2 ^y /2 <oo, 

B \\-D\ JjBx(0,r B ) 1 J 

the supremum being taken above all balls and rp denotes the radius of B. The spaces 
It-, 1 < P < 00, form a complex interpolation family. For more see |CMSj . Note that 
the IP boundedness of a conical square function reformulates canonically as an LP to 
boundedness. 

We first prove boundedness and sharpness for each square function. We consider 
next the converse inequalities globally. 

3.1. Proof of Theorem 13.11 for Gh- This was treated in [Aus] . There the range of 
p is shown to be the largest possible open set. 

3.2. Proof of Theorem 13.11 for Gh- For p < 2, it is in |Aus] . For p = 00, we first 
obtain the boundedness of / 1— > tVe~* L f from L°° to T£° by a well-known argument 
of Fefferman-Stein |FSj . More precisely, we fix a ball B and write / = f\ oc + / g i Q b 
where /i oc = f Xab- Using the L 2 boundedness of Gh, 
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Next, the off-diagonal decay (13. f I) with p = q = 2 for t 1 / 2 Ve tL implies for some 
< c, C < oo, 



— / |tVe^ 2i / glob (x)| 2 ^<L7Ve"— / 

\ t >\ JB 



\f(x)\ 2 dx 



B 



which, integrated against dt/t in t e (0,re), yields a bound by H/H^,. 

Then interpolate this estimate with the boundedness from L 2 to T 2 2 , to get bound- 
edness from LP to Tf, which is the same as the LP boundedness of Qh by rescaling 
t 2 i — y t in the integrals. 

Note that compared to [Ausj . the upper bound improves from q+(L) to oo and is 
of course optimal. As for the lower bounds, we have HG^/Hp < when p < 2. 

Hence the fact that p~(L) is optimal for Gh (see [Ausj ) implies the same for Q^. 

3.3. Proof of Theorem 13.11 for Gp. We begin with removing the V part in Gp 
when q~(L) < p < q + (L). We know that VL -1 / 2 is bounded on LP for p in this range 
|Ausj . So by vector-valued (in the Hilbert space H = L 2 {R + ,dt/t)) extension ( jGral 
Proposition 4.5.9]), we have that the 

\\G P f\\ p <C\\gpf\\ p , q.(L)<p<q + (L) 

with 

9P (f)(x)=^\tL^ Ll/2 f(x)\ 2 ^ 
Next, |HM[ Lemma 7.2] using the subordination formula 



1/2 



(3.2) e~ tL/ f = C e - — fds, 

Jo V s 

proves the pointwise inequality g P < Ccjh with 

(/■oo t.\ 1/2 

and the latter is bounded on LP for p~(L) < p < p+(L). This can be proved by 
adapting line by line [Aus} Theorem 6.1]. This also follows from Le Merdy's theorem 
[LeMl Theorem 3]. 

We conclude by noticing that p~(L) = q^(L),q + (L) < p+(L) (when q+(L) < oo). 
This finishes the proof. 

That the bounds p~(L) and q+{L) are sharp follows by the same argument as for 
Step 7 in [Avls\ Theorem 6.1]. 

3.4. Proof of Theorem 13.11 for Qp. More generally, we shall discuss here the proof 
of Proposition 11.41 for which Theorem 13.11 (4) represents the case m — 0. We shall 
treat the case m > explicitly only when p > 2, as the cases m = and m > may 
be treated by the same argument when p < 2. 

This part of the proof of Theorem 13.11 (and more generally, the proof of Proposition 
11.41) is the most involved as it does not follow from other known arguments in a simple 
way. For p = 2, this is classical integration by parts. We then present arguments for 
P-(L) < p < 2 and 2 < p < p + (L)*. That P-(L) is sharp follows by the same 
argument as for Qh- 
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3.4.1. p < 2. We present two proofs. The first one uses the recent theory of Hardy 
spaces adapted to L from |HM] or from |BZ] and the second one adapts arguments in 
|Aus] to prove instead weak-type bounds. We remark as above that the same proofs 
apply to prove (jl.4p . We shall omit details and stick to m — 0. 



Proof 1. Consider the Hardy spaces H P L defined in |HMj for p = 1 and |HMMcj for p > 
1. The H\ — > L 1 boundedness of Qp is exactly |HM} Theorem 5.6]. Then interpolation 
(see |HMMc| Lemma 4.24]) with the p = 2 case, shows H P L — > L p boundedness of Qp. 
Finally, identification of H P L with IP if and only if p~{L) < p < p+{L) proved in 
[HMMcl Proposit ion 9.1] concludes the argument. 

We mention that one can also use the abstract Hardy spaces developed by Bernicot 
and Zhao in |BZj and the interpolation further developed in |Belj . Namely it suffices 
to prove an L l estimate on some abstract atoms (that is an H Fato to L l estimate with 
H Fato as in Section 3.3 of |BZ] ) and then interpolate. By checking details and values 
(left to readers) from the clear presentation in |Be2j . one exactly finds the range for 
LP boundedness when p < 2. This theory, compared to the |HM] theory, has the 
advantage of not caring much about the "right" definition of the Hardy spaces as this 
is not needed for the purpose of interpolation. □ 



Proof 2. We proceed as in |Aus| p. 61]. We need to adapt the proof of |Aus| Theorem 
1.1] to the present situation. We take A r = I — (I — e~ r L ) N with A^ > 1 an integer 
to be chosen and follow the proof of that result with T = Qp and p~{L) < p < 2. As 
Qp is bounded on L 2 and A r satisfies off-diagonal estimates in the range (p-(L), 2] it 
suffices to show that 



(3.3) 



xER n \UdQ l :Qp(j2 h ^ X ^ >C */ 3 }| - ^ / \f( x ^ dx 



where hi = (I — A n )bi and is the sidelength of the cube Qi given by the Calderon- 
Zygmund lemma \Aus\ Lemma 1.3]. We use Chebichev and Fubini 



/ < 



9 

a 2 
9 

9 

a 2 



1 

+ — 

cr 



l \Ui4Qi 



Qp[^2,h^(x) 2 dx 



p ? 



-tL 1 ' 2 



hi{y)\ 



\Ui4Qi J J\x-y\<t 
i 

// +1 (EWy)l*v w 

[ [ (X)xR«\2 0i (2/) ItVyje't^hiiy)] 



dy dt 



dx 



t n+l 

\B(y,t)\U t 4Q l \dydt 



-tL 1 / 2 



hi(y)\ 



t n t 

2 \B(y,t)\U i AQ l \dydt 
t n t 

2 \B(y,t)\U i AQ i \dydt 



t" 



t 



Q- 



glob J 
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We estimate I\ oc . Notice that since y G 2 Qi we have that B(y, t) C 4 Qi for t < cr«. 
Then, using that the collection {2Q{\i has bounded overlapping we obtain 

4,< / f (£x^)l*v 9 , ie -"''%M|) 2 ^ 

JR" Jen v ■ / £ 

sE/ ri*v„, t e-»'" ft(W r^ 

• J2Q i Jcr i 1 

E(f7 l*V»e-""% (!/ )|^ 

v Jen JR n l 



< 

no 



< 



+ T f \tL^e- tLl/2 h t (y)\ 2 — 

J cr t JR n t 

£ r f \tL^e- tLl/ \(y)\^, 

■ Jcvi Jr" t 



where we have used the solution of the Kato conjecture |AHLMcT] to replace V y by 
L 1 ! 2 . Next we use the subordination formula (13. 2p . Minkowski's inequality and the 
change of variable t H- 1' := t 2 /As, 

\t L^e-^ /2 hi (y)\ 2 ^y /2 



oo 



c r, 

OO , POO 



< />(/■ / | (i L)VV",,,fo)|^p, 

JO V Jcrf/s JR n 1 ' 

Next we take a = - — | and use the square function estimate of Mcintosh- Yagi based 
on (tL)( a+1 )/ 2 e- iL - 

r f \tLV 2 e-t Ll/2 hi (y)\ 2 ^ tV/2 

Jcri JR n t 



< fe- s ( f / \(tL)^e- tL L- a hM 2 t- a ^) 1/2 ds 

JO ^Jcr?/sJR n J / 

\ {tL) (a + l)/2 e -t L{L -a/2 hmi 2^ d \ ^ ^ 



x </R™</0 



< r- a \\L- a /%\\ 2 < r~ a WhX = r; a ||(J - A n )b, 



where we have used jAusl Proposition 5.3] in the last inequality. To conclude we use 
that / — A r . = (I — e~ ri L ) N is uniformly bounded on L p and the CZ lemma 

i i i 

<a 2 ^ / \f(x)\*dx. 
Jr™ 
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Next we estimate / g i b- We write Cj(Qi) = 2 J+1 Qi \ 2 J Qi, j > 1. By duality we 
can take a function < H e L 2 0R™ +1 , with norm 1 such that 



' t 



dy dt 



C < ( // R „ +1 (S>**to) \tV^'\(yf dydt ^ /2 

POO /» 

<J2E 2jn \Qi\ r/ \tV y>t e- tLl/2 hMH(y,t)^ 



< 



oo „ c 

EI>"i<w(/ 

i .7 = 1 ^ 



, „ + ri/2 , . . dydt\ V 2 



t 



H(y,t) 

J2i +1 ®i 1 



2 dydt\V 2 



i j=l 

where H(y) = J °° H(y,t) 2 dt/t. We estimate 1^ by the subordination formula, 
Minkowski's inequality and the change of variable 1 1— >■ t' := t 2 /4s, 



h 3 < 



^(Q;) 



'As 



dydt\V 



JCj(Qi) 



'As 



L^e-^Kiy) 



2 dydt\V* 



< 



\VtV y e- tL hi(y) 



-'CjCQi) 



2 dydt\V 



t 



ds 



ds 



< 



G L ((I-e 



-/CjCQi) 

r? L\Nu \f„.\2 



b i ){yfdy + 



2 dydt\V 2 



t 



C ] {Q l ) 



ds 

-r? L\Ni \/\2 



^((/-e-^T^d/rdy 



C 3 -(Oi) 
1 



< 2 -i»/ 24 -^ 2 -i(5-f) / i^^r 

where in the next-to-last estimate we have used [Ausl pp. 55, 56] and the notation 
there for Gx, (the first one is here the same as Gh) and in the last one the Calderon- 
Zygmund lemma. Choosing N such that 2Af+^-n>0we obtain by Kolomogorv's 
lemma 



J 1/2 < a V 2- j 

J glob ~ " / , Z 



J2\Qi\ ess MMH{yY 



/2 



J'=l 



<a / MH{yfl 2 dy 

>UiQi 
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1/2 

dtdyx 1 / 2 



< «| U ^r /2 ( / H(y)dy 
v JR« Jo 



t 

1/2 



v JR n ' 

Gathering the estimates we have obtained for I\ oc and J g i b we conclude as desired 



/<«- 2 (/ loc + / glob )<a- 2 a 2 - p / \f( x )\ p dx = — D [ 



\f(x)\ p dx. 

□ 



3.4.2. p > 2. We shall prove a more general result, namely ( II. 4p . Let m be a non- 
negative integer and set 



Q m Am = ([[ \t^Mt 2 L) m e~ tL nf(y)\ 2 ^) 

\J J\x-y\<t 1 / 



1/2 



We begin with a series of results that are concerned with functions of L in tent 
spaces. Then, we shall deal with Q m ,p- 

Consider the notation of \Aus\ p. 10]. Let be holomorphic in /i G (cj,7t/2), 
with |0(C)| < C(l+ |C|)~ S for some s > 0, C < oo and all C G Consider for a G C, 
with > 0, ip a (0 — n+c) a $(0- Remark that 



(i+C) 

/■a 

f 1 + £-l\-a 



(i + C) a 

and since ( G £ M implies C" 1 G and arg(l + C^ 1 ) G (— /i, /i), we have that 



sup 



< e 



Consider the linear operator, a priori defined for L 2 functions and valued in T|, 

TJ = & a (t 2 L)f) t>0 . 

In the statements below, constants C are allowed to depend on the real part of a but 
not on its imaginary part. 

Lemma 3.2. For 9fta > 0, T a maps L p n L 2 to when 2 < p < p+{L) with norm 
controlled by Ce^ 3 "' for any /i G (w,7r/2). 

Proof. It is enough to consider the boundedness of T Q for the vertical norm which 
dominates the conical one, see Proposition 12.11 In this case, this follows from the 
bounded holomorphic functional calculus on LP for 2 < p < p+(L) combined with Le 
Merdy's theorem |LeM| Theorem 3]. □ 

Lemma 3.3. For 9ftct > 2 p+(L) > ^ a ma P s L p to Tf when 2 < p < oo with norm 
controlled by Ce^ 3 "' for any ji G (tu,ir/2). 
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Proof. For fixed a it is enough to consider the case p = oo as one can then complex 
interpolate from jCMSj between T 2 and T 2 °°. We claim that for any 2 < q < p + (L), 
and any ball B, 



1/2 



(3-4) (± [[ \^(t 2 L)f(x)\ 2 ^) 

\\&\ JJBx{0,r B ) 1 J 

< Ce^ Y^2-^ ma ~ n /^ \f{x)\ 2 dx^j 

We postpone the proof of the claim until the end of this subsection. Now the right 
hand side is dominated by the L°° norm of / by using 5ia > 2 p+(L) anc ^ choosing 
q < p+{L) appropriately. Then the supremum over all B of the left hand side is 
precisely the T 2 °° norm of T a f. □ 

Lemma 3.4. For < < T a maps LP to 7f when 2<p< n Jg+ jRa . 

Proof. By a result of Harboure, Torrea, Viviani |HTV] , there is a linear map i which 
for all 1 < p < oo is an isometry from Tf to a closed subspace of L P H where H = 
L 2 (R" +1 ,^f). Thus, the maps toT a form an analytic family of linear operators 
and they are bounded from LP to L P H for (l/p,a) given by the two above lemmas. 
Stein's complex interpolation theorem (see |Gra[ Theorem 1.3.7]), extended to H- 
valued functions (use the linear C-valued maps / t— > (toT a (f), h) for any fixed h G H), 
applies since the growth is controlled in and gives the desired range of p in terms 
of 9ta. □ 

We can use the above combined with the following lemma whose proof is postponed 
to Section @] 

Lemma 3.5. Let m be a non-negative number. For C depending only on ellipticity 
and dimension, for any function f e L 2 and any x € M. n , 

1 /2 

g^(f)(x) <mC(ff WLTe-« L f){y)\ 2d ^\ 

\JJ\x-y\<2t 1 J 

1/2 



c([[ KA(^re-'' L fW)\ 2d ^) 

\JJ\x-y\<2t 1 J 



+ C( I \((t 2 Lr(e-^f(y)-e^ L f)(y) 

J\x-y\<2t 



1/2 



We now conclude for Qp = Qo t p. Start from the decomposition in the previous 
lemma and notice that the first term vanishes since m — 0. The second term is 
bounded on LP for 2 < p < oo using and rescaling t i— >■ t 1 / 2 for the part with V y 
and the same argument applies for the d t part because it picks up one more power of 
L and one still has good decay in the L 2 — L 2 off-diagonal estimates. 

For the term with e~* Ll/2 — e~* 2i , we apply the third lemma concerned with T a with 
0(() = (i + £)i/2£-i/2( e -cV»_ e -C) an d « = 1/2, which gives 2 < p < n "_f+ + ( g = p+(L)* 

if p+{L) <nor2<p<ooif p+(L) > n. 

For m positive integer then the third term of the decomposition is estimated as 
above with 0(C) = (1 + C) m+1 ^ 2 C~ 1 / 2 (e _?1/2 — e~^) and a = m + 1/2, which gives 



VERTICAL VERSUS CONICAL SQUARE FUNCTIONS 17 

2 < P < n-(2m+S§ + W if (2m + < n or 2 < p < oo if (2m + l)p+(L) > n. The 

first term is as good as the second one, i.e., bounded on L p for 2 < p < oo. □ 

Remark 3.6. It seems that the order ( a for ip a at governs the p range for boundedness 
of the conical square function. But if the decay of the off-diagonal estimate is fast 
enough, then this information is not necessary. For example, consider the conical 
square function made after (t 2 L) m e~ t L for m a positive real number. When m is an 
integer, they are bounded on LP for all 2 < p < oo because the decay is gaussian 
(polynomial of some high enough degree would suffice). But when m is a non-integer, 
then the decay is polynomial and our method gives a limited range of p for small m 
unless p+(L) = oo. In other words, when p+(L) < oo, we obtain a range of p that is 
discontinuous a function of m. We do not know whether this discontinuity is a reality 
or an artifact of our method. We ask therefore whether Q m ,p is bounded on LP for 
2 < p < oo and all real m > 0. 

Proof of ([331). We write / = /loc + /glob where / ioc = f Xab- Then, using the L 2 
boundedness of square functions associated with (p a (t 2 L), 

^ [[ M 2 L)fUx)\ 2 ^<T^ [ ( r\v a (t 2 L)f loc (x)\ 2 ^)dx 

\ t} \ JJBx(0,r B ) 1 \ ti \ J9. n K J0 tJ 



It is then enough to show 



/ \fUx)\ 2 dx = C-f \f(x)\ 2 dx. 
I-d| Jin Jab 



\ 1/2 f 2dta 00 / r \ 1/2 

\ip a (t 2 L)f gloh (x)\ 2 dx) < Ce^ J2^ ji2Sta ' n/q) [} v+lB \f^\ 2dx 

J— 2 

Indeed, plugging this estimate in the integral on the Carleson region, we obtain the 
claim. 

To this end, we set fj = f Xc 3 (b) with Cj(B) = 2 ]+1 B \ 2 3 B so that / glob = Y^j>2 fj 
and by Minkowski's and Holder's inequalities 

1/2 / f \ l/q 

2rW U|2j„l / / I,. f+2r- 



^ V a {t 2 L)f gXoh {x)\ 2 dxj < ^ ^JfL^xWdx ) 



for any q > 2. Fix j > 2 and use the representations [Aus| (2.6)-(2.7)] to estimate 
<p a (t 2 L)fj. For the 77±,t(.z) given by |Aus| (2.7)] we find with v G (co,fj), 

Cf* a 
ma+i ' 



Next, using (J3HJ) in jAUsl (2.6)] for e"^ L with p = 2 and 2 < q < p+{L), E = Cj{B) 
and F = B, we easily obtain 

/ /■ \ 1/9 +25Ra / /■ \ 1/2 

(/ 1^(^)^(1)1^) < C< e ^l 2"^ 2 ^-^) -/ |/(x)| 2 rfx . 

\Jb J r B \J2i+ 1 B J 

We see in this last estimate the combined roles of and p+(L): 9ftct > yields 
integrability in t while 2§ta — n/p + (L) > yields the summability in space. □ 
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3.4.3. Converse inequalities. We basically follow |Aus| Theorem 6.1, Step 8]. What 
we have proved so far applies to any operator L in our class, and in particular, to 
L = —A. The explicit formula for the heat kernel implies that A) = 1 and 
j9+(— A) = A) = oo. Hence, we obtain the well-known estimates 

\\G P ,- A f\\ P + IIGVa/IIp + ||&.-a/||p + ||0k,-A/|| P < \\f\\p 

for all 1 < p < oo and / G U\ where we have adapted the notation to indicate the 
operator. 

The converse < is based on the following formula for /, g E L 2 : 

/ f( x ) 9~( x ) dx = lim / e~ eL f(x) e eA g(x) dx — lim / e~ RL f(x) e RA g(x) dx 



d 

dt 1 



-tL 



f){x) (e tA g)(x)dxdt 



(A(x) + /)(Ve^7)(i) • (Ve tA g)(x) dxdt. 

n x(0,oo) 

The last equality is obtained by integration by parts in the x variable after computing 
the time derivative. Hence, we obtain with obvious notation 



f(x)g(x) dx 



< (PIU + 1) / G h , L (f)G h ,„ A (g), 
it™ 



so that 



and it follows 



f(x) g(x) dx 



< ||G ! /, i L(/)||p||p||p' 



<\\G h Af)\\ 



V 



For L the proof is similar. Starting from the equality above, we use the averaging 
trick of the Introduction and then Holder's inequality. Details are left to the reader. 

For square functions based on the Poisson semigroup, the idea is the same but one 
needs to integrate by parts in t twice: 



f(x)g(x)dx = - 



00 d 
o dt 
00 d 2 



-tL 1 / 2 



f){x) (e- t (- A y /2 g)(x) dxdt 



t 



dt 2 



-tL 1 / 2 



l x(0,oo) 



f)(x) (e-*(- A ) 1/2 #)(x) dxdt 
(A(x) + I)(tV x e- tLl/2 f)(x) ■ (tV x e-^- A ) 1/2 g){x) 



dxdt 



{tV t e- tLl/2 f){x)-{tV t e^~ A ) 1/2 g){x) 



'x(0,oo) 



dxdt 
t 



The last line is obtained by distributing the second derivatives in t and integrating by 
parts in x using Jp-(e _tLl/2 f)(x) = L(e~ tLl/2 f)(x) and similarly with —A. The two 
right hand terms are controlled by both ||Gp i L/|| p ||G'p_Afi , ||p' and ||^p,L/||p||^p,-Afi , ||p' 
so that the conclusion follows as above. 
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4. Proof of Lemma 13.51 

If m = we take / G L 2 and set fo — f- If m > 1, as the domain of L m is 
dense in L 2 (R n ), it suffices to assume / in that space and we set f m = L m f. Define 
u m = L m e~ tLl/2 f = e~ tLl/2 f m , v m = L m e~ t2L f = e~ t2L f m . Notice that 

tV yjt (t 2m u m ) = 2mt 2m v m e + 2mt 2m (u m - v m )e + t 2m {tV y , t u rn ) 

with e = (0, . . . , 0, 1). The first and second terms give rise respectively to the first 
and third terms on the right hand side of the desired inequality. Therefore it suffices 
to control the third term which gives a square function that is pointwise smaller than 
the integral 

/(*)= ff\v„,u m (y,t)\V fx - At,mivii 



y,t^rn\y, "J\r I ^ I ^ n _ x j 

where ip is a smooth positive function with </? — 1 on the unit ball -6(0, 1), supported 
in the ball 5(0,2). To justify the calculations, for < r < R/10 < oo, let 4>r,n(t) = 
C(t/r)(l — ((t/R)) where £ is a smooth function that satisfies < ( < 1, ((t) = if 
t < 1/2 and = 1 if t > 2 and set 

= // |V^ m (y,t)|V(^)<«W^^. 

Let i? be the (n + 1) x (n + 1) block matrix with A being one block and 1 the other 
one. By ellipticity I r ,ii(x) < C(A)9fJX ? . jR (x) with 

2 / x — y \ 2 , N m dydt 



m— 1 



Next, we write 

'x — y\ 2 t 4m dydt 



4-n—l 



+ I/ 5(y)V^ m -V^ m ,^ 2 (^^^(t)^^=X r 1 ifl (x)+X r %(x). 

In the last integral, distribute the product on each gradient term and use Young's 
inequality with e to obtain a bound 

,t 4m dydt 



Ml— 1 



\\B\\ 00 £l rtR {x) + Ce 1 // |V y ,iV. 
Using that 

* 2m (tV y ,^ m ) = tV ytt {t 2m v m ) - 2mt 2m v m e 

we can obtain 
X 2 R (x) < H^llooe/r^a;) 

JJ|a:-j/|<2t 1 JJ|a;-j/|<2t 1 

Note that the first term can be hidden if e is small enough independently of r, R, x. 
For X} R (x) we integrate by parts using the equation satisfied by u m to obtain 

x lA x ) = ~ Jj B(y)V y , t u m ■ Vj /it |^ rT (^ 2 ^^y^^ 2 i? (t)|(n m - v m )dydt. 
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Note that 

t 4m 2 fx-y\ 2 \ t 2m (x-y\ 9(y,t)t 2m 



where 9: — > R n+1 is a function with support in the cone defined by |ar — j/| < 
It and is bounded independently of x, r, R. Hence, another application of Young's 
inequality with e yields a bound 

,t im dydt 



Tl R {x) < WBW^el^ + Ce- 1 [[ 

JJ\x-y\ 



t 



n+l 



y|<2* 

Again, the first term can be hidden if e is small enough independently of r, R, x. 
Gathering the obtained estimates we conclude that 



I r , R (x) <C(n,\,A)(m [[ \t 2m v m \ 2 ^- + [[ \V y , t (t 2m v m )\ 

V JJ\x-y\<2t t J J\x-y\<2t 



2 dydt 

t n ~ l 



JJ\x-y\<2t 1 

Letting r \. and R f oo, one obtains the desired estimate. 
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